Big Bang Nucleosynthesis with a New Neutron Lifetime 
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We show that the predicted primordial helium production 
is significantly reduced when new measurements of the neu- 
tron lifetime and the implied enhancement in the weak reac- 
tion rates are included in big-bang nucleosynthesis. There- 
fore, even if a narrow uncertainty in the observed helium 
abundance is adopted, this brings the constraint on the 
baryon-to-photon ratio from BBN and the observed helium 
into better accord with the independent determination of the 
baryon content deduced from the WMAP spectrum of power 
fiuctuations in the cosmic microwave background, and mea- 
surements of primordial deuterium in narrow-line quasar ab- 
sorption systems at high redshift. 



I. INTRODUCTION 

Big-bang nucleosynthesis (BBN) plays a crucial role 
in constraining cosmological models. It is essentially the 
only probe of physics in the early universe during the 
interval from ^ 1 — 10* sec in the radiation dominated 
epoch. As such, it is important to have accurate predic- 
tions of the light element abundances produced in this 
era. 

The single unknown parameter for standard BBN is 
the baryon-to-photon ratio during the nucleosynthesis 
epoch. All light abundances are a simple function of 
this parameter. In this regard, it has been noted for 
some time [1-7] that the nucleosynthesis yields from the 
big bang are particularly sensitive to the neutron lifetime 
which affects BBN in two ways. For one, changing the 
neutron lifetime r„ implies different weak reaction rates 
through the relation between the neutron lifetime and 
the weak coupling constant. 

^n' = 0(l + 35lKAo , (1) 

where Gi? is the Fermi coupling constant and gA is the 
axial- vector coupling of the nucleon. The quantity rrie is 
the electron mass and Aq is the phase-space integral for 
neutron decay. To a good approximation, weak reactions 
cease once the weak reaction rate 

r = (7/60)4H- 35^)G|t5 , (2) 

becomes smaller than the Hubble expansion rate. 



77«[(8/3)7rGp,]i/2 , (3) 

where — (7r^/30)(7*T'* is the energy density in rela- 
tivistic particles, and 17* the total number of effectively 
massless degrees of freedom at the relevant epoch. 

Equating these two rates gives the freezeout temper- 
ature, Tf K, 1 MeV. Once weak reactions freeze out the 
ratio of the number of neutrons/protons remains fixed at 
the freezeout value except for neutron decay. However, 
changing the neutron half life resets the temperature Tf 
at which weak-reactions freeze out. 

For example, a shorter lifetime for neutron decay 
means that the reaction rates remain greater than the 
Hubble expansion rate until a lower freezeout tempera- 
ture. This shifts the equilibrium neutron-to-proton ra- 
tio at freezeout. To a good approximation this n/p ra- 
tio is just given by thermal equilibrium to be n/p — 
exp{— Am/T/}, where Am is the mass difference be- 
tween the neutron and the proton. Since most of the 
neutrons remaining until the nucleosynthesis epoch at 
t ^ 200 sec are converted to *He, there is a simple ap- 
proximate relation between the n/p ratio at freezeout and 
the helium mass fraction from BBN 

Yp « 2n/(n -f p) = 2{n/p)/{n/p + 1) , (4) 

where n and p refer to the number densities of neutrons 
and protons, respectively. The other dependence oiYp on 
the neutron lifetime simply comes from the fact that some 
neutrons can decay in the interval between weak freeze- 
out [t ^ 1 sec) and nucleosynthesis [t ~ 200 sec). Taken 
together, both of these effects imply that the shorter the 
neutron lifetime, the lower the predicted BBN helium 
abundance. 

Regarding the neutron lifetime, it is of particular inter- 
est that a new and very accurate measurement of the neu- 
tron lifetime has recently been reported [10] using ultra- 
cold neutrons in a gravitational trap. There are several 
distinguishing features of this measurement. Among the 
most important are: 1) that it involves the best-measured 
storage time (872 ± 1.5 sec) of neutrons in the trap; 2) 
the ability to measure the spectrum of the ultracold neu- 
trons after they have been stored in the trap; and 3) an 
improvement in the coating of the traps which improves 
the reliability for the different geometries. These fea- 
tures are particularly important improvements because 
all measurements of the neutron lifetime involve an am- 
biguity between the neutron decay lifetime t„ and the 
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storage lifetime Tsurrage in the trap. To circumvent this, 
it is necessary to conduct a number of measurements with 
different storage lifetimes which can be then extrapolated 
to zero storage loss rate (Tj^oragp 0) to determine the 
decay rate due to neutron decay alone. This is a source 
of considerable systematic error. The present result is a 
major improvement in previous extrapolations in that the 
neutron storage loss rate was not only accurately mea- 
sured but was as much as a factor of two smaller than the 
best previous measurements, making the inferred neutron 
decay lifetime much less subject to systematic error. In- 
deed, the difference between the best-measured storage 
time and the inferred neutron lifetime is only 5 s, whereas 
in the previous best measurements the extrapolation was 
made over an interval of 105 s and therefore less reliable. 

The neutron lifetime deduced by this method is signif- 
icantly reduced to 878.5 ± O.lstat ± O.Ssys sec. This value 
differs from the previous mean weighted world average of 
885.7 ± 0.8 [9] by six standard deviations. It differs from 
the previous most precise result of 885.4 ±0.9sfaf±0.3sys 
[11] by four standard deviations. Indeed, including the 
new result of Ref. [10] into deriving a new weighted mean 
world average according to the methods of the Parti- 
cle Data Group [9] reduces the mean weighted world av- 
erage by over four standard deviations to 881.9 ± 0.6. 
More conservatively, however, this weighted mean uncer- 
tainty may not be appropriate due to the inconsistency 
among the data. A chi-squared minimization instead of a 
weighted mean gives a larger uncertainty of ±1.6 s, which 
we adopt here. This larger uncertainty and smaller life- 
time together help to bring the primordial helium abun- 
dance into concordance with other determinations of the 
baryon-to-photon ratio as described below. 

On the other hand, there are reasons to consider the 
new value by itself, independently of the previous mea- 
surements. In addition to the vast improvement in the 
present measiircmcnt, another aspect which lends par- 
ticular credibility to this new result is the fact that 
when this new lifetime is used as a unitarity test of 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix [10], to- 
gether with the current value of the /j-asymmetry in neu- 
tron decay [12], there is excellent agreement with the 
standard-model prcxiictions. Such is not the case for the 
current world average [10]. 

If this new lower value for the neutron lifetime is 
adopted as a most extreme case, then this substantially 
reduces the expected "^He abundance from primordial nu- 
cleosynthesis. This is particularly important for BBN 
cosmology as we now explain. 

II. LIGHT ELEMENT ABUNDANCES 

One of the powers of BBN is that all of the light ele- 
ment abundances are determined in terms of a single pa- 
rameter ?7io which is the baryon-to-photon ratio in units 



of lO"""^*^. The crucial test of the standard BBN is, there- 
fore, whether a single value of rjio can be found which 
reproduces all of the observed primordial abundances. 
The different light element abundances are determined 
by different means. This makes each determination an 
important independent check on BBN. 

Primordial deuterium is best determined from its ab- 
sorption line in high redshift Lyman a clouds. The aver- 
age of measurements of six absorption-line systems to- 
wards five QSOs gives [13] D/H = 2.78ta,tt x lO"'"^- 
This would imply an value of 7710 = 5.9 ± 0.5. This is 
an important result because it is also very close to the 
value Qbh"^ = 0.0224 ± 0.0009 (7710 = 6.13 ± 0.25) de- 
duced [14] from the WMAP independent determination 
of the baryon content at the epoch of photon last scatter- 
ing. Because of the concordance of these two independent 
methods, the WMAP determination of 7710 is generally 
accepted as the most accurate determination. 

The primordial lithium abundance , on the other hand, 
is inferred from old low-metallicity halo stars. Such stars 
exhibit an approximately constant ("Spite plateau") 
lithium abundance as a function of surface temperature. 
This is taken to be the primordial abundance. There is, 
however, some controversy [15,16] concerning the deple- 
tion of ''Li on the surface of such halo stars and/or dur- 
ing the big bang itself [17]. For the present purposes we 
adopt the value from [18] ■^Li= 1.23lgJ^ x 10"^°, where 
the errors are 95% confidence limits. 

The primordial helium abundance is obtained by mea- 
suring extragalactic HII regions in low-metallicity irreg- 
ular galaxies. Often in the past, the primordial helium 
abundance Yp so deduced tended to reside in one of two 
possible values (a low value, e.g. Yp = 0.238 ± 0.002 ± 
0.005, [20] and a high value yp0.2452±0.0015 [19]). There 
is also, however, a current dilemma regarding the uncer- 
tainty in the observationally determined primordial he- 
lium abundance. Many recent evaluations (e.g. [19]) give 
a rather narrow range of abundance uncertainty. For our 
purposes we adopt the value of [19] as a representative 
result. On the other hand, the extent of systematic errors 
in these analyses is still being debated. Another recent 
study [21] has adopted a more conservative approach and 
concluded that correlations in various uncertainties could 
stretch the error in the inferred primordial abundance. 
Their represantitvc analyses yields Yp = 0.249 ± 0.009 
and they argue in favor of range of allowed values of 
0.232 <Yp< 0.258. 

While this is being sorted out, however, it has been 
deduced by several authors (cf. [17,22]) that the com- 
bined deuterium and WMAP constraints on the baryon- 
to-photon ratio implies that the primordial helium abun- 
dance should be Yp = 0.2484g:ggg| [22] or Yp = 0.2479 ± 
0.0004 [17]. 

If we adopt the narrow helium abundance of [19] and 
the WMAP constraint of [22] there is, therefore, a possi- 
ble 2-3a discrepancy between the ^He -|- ^Li and the D 
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+ WMAP results. This dilemma with regards to BBN is 
depicted by dashed lines on Figure 1. 

One of course could (and probably should) disregard 
this dilemma if the uncertainty is as large as deduced in 
[21]. However, if this dilemma is real, then it may provide 
insight into new physics beyond the minimal BBN model, 
for example, branc- world effects [23], cosmic quintessence 
[24], time varying constants [25], etc. [4]. In this paper 
we point out an important result, however, that even 
if the most narrow uncertainty in the deduced primor- 
dial helium is adopted, then a significant portion of the 
discrepancy between BBN and the CMB results can be 
accounted for simply by adopting the new neutron life- 
time. 

III. RESULTS 

For illustration of the implications of the new value 
for the neutron lifetime, we have made calculations of 
standard homogenous big bang nucleosynthesis for three 
values of the neutron lifetime. These are: 1) the previous 
world average (885.7 ±0.8 sec); 2) the new world average 
(881.9 ± 1.6 sec) which includes the new measurement of 
Ref. [10]; and 3) the newest lower value of [10] (878.5 ± 
0.7 ±0.3 sec). 

The benchmark code used for the present illustration 
is the standard big bang nucleosynthesis code originally 
developed by Wagoner [1] and made user friendly by 
Kawano [26]. This code is available for public down- 
load [27]. The reaction rates and uncertainties are those 
adopted in [5]. Although newer reaction rate compila- 
tions and uncertainties have been evaluated [6-8], this 
code is readily available and adequate for the benchmark 
comparison of interest here. 

Figure 1 compares the primordial nucleosynthesis 
yields based upon both the previously adopted world av- 
erage with the yields based upon the new neutron life- 
time measurement and its uncertainty. The insert shows 
an expanded view of the primordial helium abundance 
for 7^10 values near those allowed by the various observa- 
tional constraints. 

From this figure it is clear that the primary effect of 
altering the neutron lifetime is to lower the primordial 
helium abundance prediction. The uncertainty in the 
predicted Yp is indicated by parallel bands on the figure. 
The uncertainty remains the nearly same with the new 
lifetime because the uncertainty in the new neutron life- 
time is nearly the same as that of the previous world av- 
erage. Though not shown on this figure, the uncertainty 
in predicted Yp increases by a factor of « 1.5 if the larger 
error (±1.6 s) in the new world average is adopted. 

The effect on other light elements is so small (^;1%) 
as to be indiscernible from the line widths on the fig- 
ure. The key point of Figure 1 is that now the pri- 
mordial helium abundance required for the baryon-to- 



photon ratio given in the WMAP and/or D/H QSO 
absorption- line results reduces from Yp = 0.2479 ±0.0006 
to 1^ = 0.2463 ± 0.0006 when using the new value for 
new neutron lifetime. For comparison, incorporating the 
new lifetime measurement into a new mean world av- 
erage would require Yp = 0.2470 ± 0.0009. These two 
later values overlap (within la) with the uncertainty of 
even the narrower of the observationally inferred helium 
abundance [19] of Yp = 0.2452 ± 0.0015. 

Alternatively, the 7710 values implied by an observed 
helium abundance of Yp = 0.2452 ± 0.0015, are 5.5 ± 0.9 
for the new lifetime, or 5.1 ± 1.1 for the new weighted 
mean lifetime as compared to 4.8 ± 0.8 based upon the 
previous world average. These are to be compared with 
the WMAP + D/H determination of r?io = 6.13 ± 0.25. 
Hence, even with this small correction to the neutron life- 
time, and adopting a narrow range for the observational 
uncertainty in Yp, the implied r/io value for either the 
new lifetime or new weighted average now overlaps the 
value required by the WMAP and D/H analysis. This 
significantly further constrains nonstandard models for 
BBN, and strengthens the viability of standard BBN as 
a probe of cosmology. 

Of course, one must still deal with the problem of ''Li 
overproduction in BBN which will have to be resolved 
by ''Li destruction, either within the big bang itself [17] 
or during subsequent stellar evolution [15,16]. We also 
emphasize that there is still additional uncertainty in 
the BBN production of helium and other light-element 
abundances due to uncertainties in nuclear reaction rates, 
particularly the d(p,7)^He, d(d,n)^He, d(d,p)^H, and 
3He(a,7)^Be. 
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FIG. 1. Predicted BBN light- element abundances vs. the 
baryon-to-photon ratio ?7io in units of 10"^*^. These are com- 
pared with the observationally inferred [17] primordial abun- 
dances (horizontal linos) and the independent determination 
of 7710 from the WMAP results (light shaded region). The 
top box shows the primordial helium abundances. The insert 
shows an expanded view of Yp near the allowed region. The 
banded regions indicate the range of predicted Yp due to the 
neutron lifetime uncertainty. The upper lines are based upon 
the previous world average r„ — 885.7 ±0.8 s. The lower lines 
are based upon the new measured value of t„ — 878.5 ± 0.8 s. 
The previous allowed 7710 values (shown by the dashed open 
box) shifts to the dark shaded box if the new neutron lifetime 
is adopted. 
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